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1 Introduction 

Recently, new analytical results pertaining to the spectrum of 2 + 1 dimensional Yang-Mills 
theory (YM2+1), that are in excellent agreement with the lattice data in the planar limit PP, 
have been derived in 0. These analytic results regarding the mass gap, string tension and 
the glueball spectrum resulted from a determination of the ground state wave-functional 
in the planar limit 3 . The analytic calculations in YM2+1 were based, as in ^, on the 
remarkable work of Karabali and Nair |S] . 

Because in 2 -t- 1 dimensions space is two dimensional, one can work in a complex basis 
in the Hamiltonian picture. One might wonder whether a similar success can be achieved in 
the more realistic case of a 3 + 1 dimensional Yang-Mills theory (YM3+1). As we will explain 
in this letter, a formalism does exist which closely parallels that of KN. The ingredients of 
this formalism were introduced long ago by L Bars |H] in his work on local gauge invariant 
'corner variables'. (Other related work can be found in [7].) 

In this rather programmatic letter, we outline and extend the general formalism, and 
explain the physical interpretation and use of corner variables both in YM2+1 and YM34-1 . In 
particular, in the context of YM2+1, the corner variables are appropriate to a real coordinate 
basis as opposed to the complex basis of KN. Also, as pointed out by Bars in |B], the 
corner variables can be also used in a covariant Lagrangian picture, which is an obvious 
additional attractive feature of the formalism. In the Hamiltonian picture the use of corner 
variables also reveals unexpected parallels between pure Yang-Mills theory in 2 + 1 and 3-1-1 
dimensions. We concentrate here on the discussion of the vacuum wave- functional in the 
Hamiltonian picture, especially in light of the large N lattice simulations. 

The themes outlined in this paper are explored in more detail in a companion paper [5] 
and in ^. The organization of this letter is as follows: in Section 2 we discuss the formalism 
of corner variables in view of our recent work on pure Yang Mills theory in 2+1 dimensions 
0. Then we turn to the more dynamically relevant issues in Section 3. Section 4 is devoted 
to the general analysis of the structure of the vacuum wave-functional. We conclude the 
paper with a short outline of some of our current work in progress. 

2 The Formalism of Corner Variables 

The key insight we start from is that the KN variables 5 in the 2 + 1 dimensional setting 
are related to line integrals in from infinity to a point x. The analogous variables in 3 + 1 
were introduced by Bars and we call them Mi{x) where Mi{x) are unitary matrices.^ They 
satisfy the defining equation 

Ai = ~d.iMiM~^ (no sum on i) (1) 

In other words 

M,{x) = Pexp[- I A] (2) 



'^We use notation similar to that of KN and warn the reader of notation difference compared to Bars. 



Figure 1: The 'corner variable' Hij{x). 



where the integral is a straight spatial contour for fixed x^ ,i ^ j. The interesting thing about 
these variables is that the natural lattice formulation parallels the continuum formulation. 
The above formulation is appropriate to the Hamiltonian formalism, if j is a spatial index. 
The translation to 2+1 KN variables is then just 



M, =M, 



M,_ 



Mt-' 



(3) 



Note that if we had used a real coordinate basis, then we would have had a pair of (unrelated) 
unitary matrices Mi, M2. Next, define the corner variables 



H,, = MrM, 



(4) 



Note that Hjj — 1 and Hji = H^, - this just means traversing the corner in the opposite 
direction is precisely the inverse group element. The Hij are unitary in a real coordinate 
basis. There is also a constraint (here written for 3 + 1 — it is trivial in 2 + 1) 



HijHjkHki — 1 



(5) 



Another description is in terms of a "semi-complex" coordinate basis {u, z, z}; then one 
may parameterize as 



H,, 



H, H, 



H\ H,, = H^H. 



(6) 



For example, one could use the notation M^ = M, Mg = M^ , Af^M„ == 1, with H = 
M^M. The constraint takes the form H^uHuzHzz = 1- In other words, there is in I? = 3 
in the semi-complex coordinate basis a complex _ff-field (compared to a Hermitian field in 
D = 2)] thus there are twice as many degrees of freedom, as expected. Returning to the 
general notation, gauge transformations act as 



Mi 



gM, 



(7) 



so the Hij are gauge invariant. 

The formalism of Bars also contains the notion of (or a generalization of) holomorphic 
invariance (which is familiar from the 2 + 1 KN setting). Indeed, note that one may introduce 
'currents' 



Jy - {djHij)H, 



-1 

ij 



(8) 



(in 2+1, J r^ Jzz and ,P ^ ~Jzz-) The extra symmetry acts as 

M,^M,h^\x^), j^i (9) 

(this apparently, as far as we can tell, was not pointed out in [5].) The condition j ^ i on 
the function hi is the analogue of holomorphy. This leaves the gauge fields invariant, and 
one finds 

Hy ^ h,H,jhj^ (10) 

In the semi-complex basis, we would have M„ — > Mu h~^{z,z), Mz -^ Mz h^{u,z) and so 
H -^ hu{z, z)Hh^{u, z), W —> h{u, z)H^h~^{z, z). The Jij transform as connections 

Jy i~^ h^J^jh^^ + djhi /i^i (11) 

In the real coordinate basis, it appears that there are six currents that are apparently 
distinct. However, there is a 'reality' condition on their derivatives of the form 

d,J,, = -H,j(d,J,^H;^ (12) 

(in D = 2, there is a similar relation which reads dJ — H{dJ^)H^^ .) By defining Jij — 
—HijJjiH~: , we may rewrite this as 

OiJij ^ CljJij — [Jij,Jij\ l^'jj 

and so there are covariant derivatives Dij = dj — Jij. These currents are related to the 
magnetic field Fij — diAj — djAi + [^i, Aj] by 

d^J^j = -Mr^FjM,. (14) 

As a short form then, we will denote 

Bi^i = diJi.i+i (cyclic). (15) 

These fields transform homogeneously under the 'holomorphic' symmetry 

S,;_i ^ h,B,^ihr^ (16) 

Note that there are only really two of these fields, because, solving the constraint ©, there 
are only two independent iJy 's, say H12 and H13. 

Note that the covariant derivatives can be written in terms of the usual derivative and 
Mi as 

V., = d, + A, = MAM~\ (17) 

Finally we note some notations and conventions used in this paper. The connection 
is expanded in terms of anti-hermitean generators T^, Ai = A'^Ta satisfying the algebra 
[Ta, Tb] — fab'^Tc- Also, we denote by tr the trace in the fundamental representation, so that 
ir(l) = N. This trace is normalized as ~-2tr{TaTb) — Sat- In the adjoint representation the 
generators are given by (T°')iic — —faba where indices are raised or lowered with the metric 
Sab and the trace in the adjoint is denoted by trad- A group element M is represented 



in the adjoint by Mat = —2tr{TaMTi,M ^), clearly we have (Af ^)ab — M^a and also 
MTbM-^ = TaM^b- 

We will need a variety of Green's functions. We denote by Gi{x,y) the inverse of di, 
such that dfGi{x,y) — 5'^'^\x^y) (no sum). Since di admits zero modes, its inverse is not 
uniquely defined and we will work with the expHcit choice Gi[x) = 9{xi)5{x2)5{x3) and 
Gi(x,y) = Gi{x — y) (which is not antisymmetric). This choice is not arbitrary, it is the 
unique choice consistent with the definition of the variable Mi as an ordered exponential. 
Indeed we can write 

M,{x) = 5](-l)" I dy{G,A,r{x,y), (18) 

where iG.iAi)^{x,y) = J dzG^{x,z)A^{z)G^{z,y)Ai{y) . 

3 Further Technical Details 

The Yang-Mills action is taken to be Sym = y^ J t^iFfiiyF^") and the Hamiltonian in the 
Yang-Mills variables is 



-i:{4(^)'-^'^")'} 



(19) 



with F° — \^ijkF% and should be supplemented by the gauss law constraint Vi^^ = 
0. The wave-functionals in the Schrodinger representation of pure Yang-Mills are gauge 
invariant functionals of A^ and the scalar product is given by 

m?^! Dii{A)nA]nAi (20) 

JA/g 

Where the integral is over the space of gauge connections modulo gauge transformations, 
D^{A) = y^wg, = voifg-)' ■ From the previous section we know that we can equivalently 
describe gauge invariant wave functionals as 'holomorphic' invariant wave-functionals of Hij 
or Jij . This change of variables can be explicitly performed both at the level of the measure 
and the Hamiltonian H; we recall here some of these results. 

Since Ai — —diAliM"^, SAi = —{'ViSAIi)AI^^, the change of variables involves a deter- 
minant 

e^ = det I '^^' I = det(ViV2V3). (21) 

The variational derivative of the action is found to be trivial 

ST 



= trad[iV,)-\x,x)Ta]=l. (22) 



SAf{x) 
The scalar product written in term of the Bars variables is simply 



*|P = J DH12DH13 nH]^[H] = JUdH^j 6{Hi2Hi3H23) *[i/]*[ff] (23) 



where DH denotes the product over left-right Haar measures on SU{N). In the last equality 
we inserted a delta function on the group to emphasize the symmetric form of the measure 
under permutation of indices. The interpretation of this delta function constraint should be 
clear: it is the integrated version of the Bianchi identity e^^'^ViFjk — expressed in terms 
of gauge invariant observables. One important feature of this scalar product is the fact that 
the identity functional '^[H] = 1, which can be can be viewed as a limit of the identity 
cylindrical functional, is a normalisable wave-functional. 

The potential term of the Yang-Mills Hamiltonian is easily expressed in terms of the 
Bars variables 

lT.(Pt)' = ^E/(^»^M+i)'(^)d2^- (24) 



^-2. 

t.a i.a 



and the kinetic term reads 



a.i 



where P^ are right derivatives on the group 

K(x), M,{y)] ^ H Ah{y)Ta 5(^\x, y) (26) 

and we have introduced the kernel 

Qf{y,z)^5'^\G,G,){y,z) = 5'''' f dx G,{y,x)G,{x,z). (27) 



It is important to note that the form (|25|) for the kinetic term is valid only when T acts on 
holomorphic invariant states. This Hamiltonian can be checked to be hermitian |S|. 

Two technical points are in order here. First, in the KN formalism, a non-trivial Jacobian 
appeared and is given by the exponent of a gauged WZW action. This is an artifact of the 
complex polarization adopted in the formalism, and it is not a universal feature. In the 
real basis that we have discussed here, such a Jacobian does not arise. Indeed, consider the 
computation of such a Jacobian determinant, which one can write conveniently as a fermionic 
determinant. Depending on the choice of a chiral (in the complex KN-like polarization) or 
non-chiral basis (in the real Bars-like polarization) for the auxiliary adjoint fermions, in the 
case of 2 + 1 Hamiltonian formulation, one gets the usual anomalous term represented by 
the WZW action, or one gets a unit Jacobian, respectively. Thus one also has a trivial unit 
Jacobian in 3 -I- 1 dimensions if one works with Bars' corner variables. 

Second, there is a delicate issue of holomorphically invariant regularization procedure 
(which extends the techniques used by [1]). And indeed, in order to make sense of the 
computation of the determinant or the hermiticity property of the kinetic term, one needs 
a regularisation scheme which preserves all the symmetries. Such a regularisation scheme is 
fully described in |H] and is absolutely crucial for actual computations. 

3.1 Semi-complex coordinates 

One of the most important observations about the Bars' formalism is that the kinetic terms 
have the form T = J{Gpq)^ (where p and q denote the canonical generalized momentum 
and coordinate, [q,p] ~ i) with a regularized Green's function G. 
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The significance of this will be appreciated by considering a toy example with a similar 
Hamiltonian. Consider a free Id rotator. The kinetic term is T = —-gffi- Here, is a periodic 
variable and the spectrum is discrete, with eigenvalues n^, and with periodic eigenfunctions 
gine ^ ]\[q^ jf -^^rg change variables p = e*^ , the kinetic term becomes homogeneous 

'^=^PTPh (28) 

op op 

where obviously p and its conjugate momentum j- have the canonical commutation rela- 
tions. So this has a (qp)^ form, as required. Now we can rewrite this using the canonical 
commutator as 

which looks like the usual kinetic term for pure Yang-Mills theory in terms of complex corner 
variables as in [Hj. We look at the action of this kinetic term on p", rt = 1,2, 3... which 
correspond to the eigenstates e*"^. The first term gives a homegeneity factor (the Euler 
factor) of n, which is of course not the correct eigenvalue! To get the correct eigenvalue n^ 
we need also the second term in T^. 

The upshot here is that we see a close parallel between this simple toy model and the 
formulation of pure Yang-Mills theory in terms of corner variables. 

In view of this toy example we discuss the complex versus real corner variables and the 
form of the corresponding Hamiltonian. The similarity between the Karabali-Nair Hamil- 
tonian written in terms of complex corner variables and the collective field theory |10| has 
been noted before [1], ^, [21 ■ Yet the Hamiltonian written in terms of real variables ap- 
parently does not have this form. In particular, it seems that the homogeneous piece of the 
Hamiltonian (as illustrated by the above toy example) is missing. One would have such a 
term however, if we write the theory in the semi-complex basis in 3 -f 1 that we discussed 
in Section 2. Such a term acts homogeneously on suitable functionals and thus acts as an 
effective dimension counting operator. In what follows we will further comment on this issue 
when we discuss the constituent picture of glueballs. 

More explicity in terms of the semi-complex coordinates from section 2 (u, z,z), the 
gauge invariant information is encoded in a complex H = Huz- In this case the variational 
derivative of the determinant ^^j", -, can be computed to give |H] 

ST = -^ ftr[{H^H)-'5{HH^)d{iH^H)-^dH^H)] (30) 

T:-i/^ J 

This expression integrates to a 3 dimensional generalization of the WZW action 

2Np - 



-3/2 



(31) 
where B^ denotes a 3-ball bounding the plane u = const. 

In this form the parallel with the KN formalism j2| becomes very striking. The 2 + 1 
dimensional results developed in j2| can be thus viewed as a natural reduction of the 3 + 1 

^Note that in the quantum field theory setting the factor of n^ is misleading as it can be seen in the 
detailed analysis of the 2 + 1 dimensional Yang-Mills theory |2]. Further discussion is in Section 4. 



dimensional formulation in the semi-complex basis. In particular, the Hamiltonian in the 
semi-complex basis is |9j (focusing on the Jzz dependence) 



T-t ~ — :rj:T { / dud z J^zTl ^ du d wid W2 ^{wi,W2)-f-r-, T-rr~^ T + • 

27r''/^ U OJzz J J OJzz[Wl) dJzz{w2) 

(32) 
where H. is essentially fixed by the two-point function of the WZW model as in the 2 + 1 
dimensional context [51 12] ■ 

4 Vacuum wave-functional: general discussion 

In this section, we make some general comments about the vacuum wave-functional in 3 -I- 1 
dimensional Yang Mills theory in view of the already discussed parallel between pure gauge 
theories in 2 + 1 and 3 + 1 dimensions. First, on general grounds one knows that the vacuum 
wave-functional 'So is gauge invariant, parity even (in the absence of 9 terms) and most 
importantly should be a strictly positive functional. This can be written as 

*o = e^ (33) 

where P is a functional of the 'position' variables. Second, we also know that the vacuum 
wave-functional can be formally written as a path integral on half space-time as 

*o[i] = / DA e^^--''^"'^ *^(F.,.f-)(M) (34) 

JA(t=0)=A 

where we integrate over all fields satisfying the boundary condition A(0, x) — A{x) and the 
integral is computed in the Aq = gauge. This implies that P is a gauge invariant sum over 
connected Feynman diagrams and in the large N limit is a non-local functional involving 
only a single trace over gauge indices. 

A physical mass scale m must emerge in this theory. In fact, in order to write a wave- 
functional that involves arbitrary number of derivatives acting on fields, a scale must be 
introduced by hand. Presumably then, physical mass scales would be determined self- 
consistently. Furthermore, if the 3 + 1 theory follows the 2 + 1 theory, we can suppose that 
in the IR, the mass scale enforces locality of P in this limit. Ultimately, this question can only 
be answered for sure by computing physical quantities (such as the Wilson loop expectation 
value). Given this assumption, we can formally develop the non-local functional P in terms 
of local operators. By the previous argument, this expression involves only local, single 
trace gauge invariant functionals. Such operators are naturally labeled by two integers n, k 
which count the number of laplacians and powers of magnetic field insertions respectively. 
We can write this expansion, in parallel with 2+1, symbolically as 



^ k,n '' 



4) (4 



(35) 



We will advocate, as in 2 + 1, that it is useful to take P to have a certain quadratic form. 
In particular, when written in terms of the gauge invariant variables, we will consider as an 



ansatz 

*o = exp (^ f tr B, K{L/m^)B, + . . .V (36) 

where X is a kernel constructed out of 'holomorphic' covariant derivatives scaled by to, 
but does not contain the curvature Bi. The kernel K should be determined (although we 
do not do so here) consistent with gauge, 'holomorphic' and spacetime symmetries. Most 
importantly, it should be determined such that the wave- functional is normalizahle. In 2 + 1, 
such a kernel was found, and it seems reasonable to suppose that this could be repeated in 
3 + 1, given the similarity of the formalism. 

We note, independent of these remarks, that P should have a specific asymptotic form in 
the far UV — the theory is free in that limit, and thus the vacuum wave-functional should 
be appropriate to free gluons. In this limit, we should find 

Puv^\ I tr B,{-k)-^B,{k). (37) 



This is determined essentially by dimensional analysis and is consistent with the transver- 
sality of gluons. We note that the mass scale m introduced to define the kernel, disappears 
in this limit. It is a straightforward exercise to demonstrate that the Schrodinger formalism 
does indeed reproduce this perturbative result. 

The UV behavior of the wave-functional is required but does not address the question of 
normalizability Also, to be consistent with confinement, (in the sense of no phase transition 
occurring between UV and IR) the kernel K should behave, thought of as a function of 
momentum, in a smooth way, and as we have discussed, take a certain simple form in the 
far IR. In particular, we are supposing here that this form is 

PiR = tA- I tr B,{-k)B^{k). (38) 

^5 ™ Jk 

Another way to motivate this form would be to compute how the Hamiltonian acts on this 
operator, and build the Schrodinger equation. Using the 'holomorphic' regularization [5] 
with a regulator scale fi and in the context of real variables, one finds 

g^T ■ j BfBf = 2M J BfBf (39) 

where 

2M=^.. (40) 

This is reminiscent of the behavior in 2+1, and is consistent with the kinetic energy acting 
to 'count' derivatives. Perhaps the simplest way to understand this result is to resort to the 
semi-complex basis in which the intuition gained in 2 + 1 dimensions becomes useful. In 
that basis, there is a homogeneous part in the kinetic term (as in our toy example) which 
acts as a dimension operator. 

H^M f J,,-^ + ... (41) 



where M is an effective mass scale. This is a first order operator. The rest of the kinetic 
term, a second order operator, acts to properly normal order the operators in the regularized 
calculation as well as provide the necessary invariant counterterms so that the final result 
is indeed holomorphic invariant. 

Note though that the difference between 2+1 and 3 + 1 is in the power-counting. In 
2 + 1, M is naturally given by the coupling constant due to its dimensionality. In 3 + 1 the 
coupling constant is dimensionless, and M is generated dynamically from the regularization 
of T, and appears for obvious dimensional reasons! It becomes clear though that this result 
is regulator dependent, and what happens in the continuum limit is not clear. 

The result H39|) is consistent with (|38|l . relating M to m. The IR vacuum wave- functional 
provides a probability measure ^q^o equivalent to the partition function of the Euclidean 
three-dimensional Yang-Mills theory with an effective Yang- Mills coupling ^g^ = mg^ . 
One could use this to compute the expectation value of a large Wilson loop and deduce the 
area law behavior. Presumably this would mean that the square root of the string tension 
scales with m. This remark ties the scale introduced into the vacuum wave-functional to 
this physical parameter. 

Let us remark further on the generalized Gaussian form of the vacuum wave-functional. 
In particular, the neglect of higher order terms in the exponent of the wave-functional needs 
to be justified. We want to be perfectly clear that the validity of the generalized Gaussian 
form for the vacuum wave-functional cannot be established by appealing only to standard 
large N simplifications. The large N limit only selects single trace expressions in the ansatz 
for the wave- functional. Thus we also do not expect that this wave-functional is exact; if 
this were true, Bi would indeed represent the right variables appropriate to the large N 
limit. Nevertheless, as in the 2-1-1 dimensional case 0, these local gauge invariant variables 
are the correct constituent degrees of freedom, even though they do not appear as physical 
asymptotic states! The physical states (an infinite number of glueballs, or non-local gauge 
invariant variables, such as Wilson loops) can be built out of these local degrees of freedom, 
so that all expected predictions based on the large N counting (factorization, suppression 
of vertices etc) are fulfilled. 

The neglect of the higher order terms in Bi suggests that there is a second expansion 
parameter at work, involving a new length scale. One such candidate scale is the size of 
glueballs. The approximation employed here amounts to considering only 'free' glueballs, 
which are point-like and non-interacting. This is consistent with the large N picture where 
one expects that all glueball interactions are suppressed by powers of 1/A^. This wave- 
functional has the form of a "generalized coherent state" appropriate to large N J^ . One 
way to intuitively motivate the Gaussian form of the wave-functional is to think of Bi as the 
relevant local probes of real physical states, but not as actual asymptotic states! Then apart 
from the rank of the gauge group N , there should exist another expansion parameter, which 
is related to the size of the glueballs. The quadratic term in the wave functional should be 
then interpreted as the leading term in the expansion in the inverse of that effective glueball 
size. This would be very reminiscent of the a' expansion in string theory. 

As we mentioned above the confinement should be implied by the normalizability of the 
wave-functional, as in the 2-1-1 dimensional counterpart T . Note that the constituent Bi 
should not appear as an asymptotic state. Also, we further expect the glueball constituents 
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to be "seeds" for constituent quarks once the fermionic degrees of freedom are included. 

Furthermore we also want to reiterate that the "QCD scale" as given by m should be 
self-consistently determined: on one side m is needed in the wave-functional for dimensional 
reasons, and thus it sets the mass scale for the spectrum, and on the other hand, the square 
root of the string tension is given in terms of ?7i up to a multiplicative numerical constant. 
Note that the cut-off /i appears both in the expression for the gap and the string tension, 
and should consistently cancel in the ratio (assuming they both persist in the continuum). 
Also, in the limit of small coupling, the ratio of the string tension and the scale m (or 
equivalently the cut-off /i) should only depend on the dimensionless coupling g'^N, and thus 
given the fact that the wave-functional describes a free theory in the UV, is then corrected 
logarithmically, as implied by asymptotic freedom. 

Finally, one of the amazing features of large N lattice simulations ^21 is that the actual 
numerical values for the ratio of masses and the square root of the string tension are of the 
same order of magnitude, both in 2 -t- 1 and 3-1-1 dimensions. The actual numbers are very 
close, up to 15% to 20%. We think that this is evidence that the large N Yang Mills theories 
in 2 -|- 1 and 3 + 1 theories are "close" in the sense of the respective continuum limits. We 
think that this is not a coincidence in view of the formalism we have discussed. 

5 Conclusion 

In this short programmatic paper we have discussed an analytic approach towards the 
solution of pure Yang-Mills theory in 3 -I- 1 dimensional spacetime. Our approach is based 
on the use of local gauge invariant corner variables in the Schrodinger representation and the 
large N limit. In particular, within this approach one finds unexpected parallels between 
pure Yang-Mills theory in 2 -I- 1 and 3 + 1 dimensions. 

There are obviously many important questions to be addressed. Perhaps the most im- 
portant is to determine a definite form of the kernel in the quadratic wave-functional in 
parallel with the successful discussion of the 2+1 dimensional situation 2 which did lead 
to the computation of the large N glueball spectra/Regge trajectories. The large N lattice 
results are already available jl2j. 

There are of course many other obvious questions, both in the 2 + 1 and 3 + 1 dimensional 
contexts: the inclusion of fermions and the computation of the meson (and baryon) masses 
in the large N limit; the development of the covariant approach; the elucidation of the role 
of topology; the proof of confinement and many others. We plan to address some of these 
questions in 0. 
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